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ABSTRACT 

We calculate as a function of orbital distance, mass, and age the theoretical spectra and orbit-averaged 
planet/star flux ratios for representative wide-separation extrasolar giant planets (EGPs) in the optical, near- 
infrared, and mid-infrared. Stellar irradiation of the planet's atmosphere and the effects of water and ammonia 
clouds are incorporated and handled in a consistent fashion. We include predictions for 12 specific known EGPs. 
In the process, we derive physical diagnostics that can inform the direct EGP detection and remote sensing pro- 
grams now being planned or proposed. Furthermore, we calculate the effects of irradiation on the spectra of a 
representative companion brown dwarf as a function of orbital distance. 

Subject headings: planetary systems — binaries: general — ^planets and satellites: general — stars: low-mass, 
brown dwarfs — radiative transfer — molecular processes — infrared: stars 
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1. INTRODUCTION 

To date, more than 110 EGPs (Extrasolar Giant Planets) have 
been discovered by the radial-velocity technique around stars 
with spectral types from M4 to Fl*. A sample of relevant refer- 
ences to the discovery literature, by no means exhaustive, in- 
cludes Mayor and Queloz (1995), Marcy and Butler (1996), 
Butler et al. (1997,1999), Marcy et al. (1998,1999), Marcy, 
Cochran, and Mayor (2000), Queloz et al. (2000), Santos et 
al. (2000), and Konacki et al. (2003). These planets have min- 
imum masses (irip sin(/), where / is the orbital inclination) be- 
tween ~0.12 Mjand ^15 Mj(Mj= one Jupiter mass), orbital 
semi-major axes from ^0.0225 AU to ^5.9 AU, and eccentrici- 
ties from ~0 to above 0.7. Given our all-too-narrow experience 
within the solar system, such variety and breadth was wholly 
unanticipated. 

Importantly, two EGPs (HD 209458b and OGLE-TR-56b) 
have been found to transit their primaries (Henry et al. 2000; 
Charbonneau et al. 2000,2001; Brown et al.2001; Konacki 
et al. 2003; Torres et al. 2003). Furthermore, in the first mea- 
surement of the composition of an extrasolar planet of any kind, 
Charbonneau et al. (2002) detected sodium (Na-D) in their HD 
209458b transit spectrum. This was followed by the detection 
of atomic hydrogen at Lyman-a and the discovery of a plan- 
etary wind (Vidal-Madjar et al. 2003; Burrows and Lunine 
1995). With both transit and radial-velocity data, an EGP's 
mass and radius can be determined, enabling its physical and 
structural study. Such data can resolve the ambiguity inherent 
in the radial-velocity technique's sensitivity to only the combi- 
nation nip sin(i). Precision astrometry can also be used to derive 
masses, as has been done for GJ 876b with the Fine Guidance 
Sensors on HST (Bennett et al. 2002), and space interferom- 
etry using SIM (Unwin and Shao 2000) promises to provide 
unprecedented astrometric masses by the year ^2010. 

However, it is only by direct detection of a planet's light us- 
ing photometry, spectrophotometry, or spectroscopy that the de- 
tailed and rigorous study of its physical attributes can be con- 
ducted. By this means, the composition, gravity, radius, and 
mass of the giant might be derived and the general theory of 
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EGP properties and evolution might be tested (Burrows et al. 
1995,1997; Marley et al. 1999; Sudarsky, Burrows and Pinto 
2000; Sudarsky, Burrows, and Hubeny 2003 (SBH); Baraffe et 
al. 2003). For close-in EGPs, we can anticipate in the next few 
years wide-band precision photometry from MOST (Matthews 
et al. 2001), Kepler (Koch et al. 1998), Corot (Antoneflo 
and Ruiz 2002), or MONS (Christensen-Dalsgaard 2000) that 
will provide the variations of the summed light of the planet 
and star due to changes in the planetary phase. In the mid- to 
far-infrared, the Spitzer Space Telescope (a.k.a. SIRTF, Space 
InfraRed Telescope Facility; Werner and Fanson 1995) might 
soon be able to measure the variations in the planet/star flux ra- 
tios of close-in EGPs (SBH). For wide-separation EGPs, it is 
necessary to measure the planet's light from under the glare of 
the primary star at very high star-to-planet contrast ratios. To 
achieve this from the ground, telescopes such as the VLT inter- 
ferometer (Paresce 2001), the Keck interferometer (van Belle 
and Vasisht 1998; Akeson and Swain 2000; Akeson, Swain, 
and Colavita 2001), and the LBT nulling interferometer (Hinz 
2001) will be enlisted. From space, the Terrestrial Planet Finder 
(TPF, Levine et al. 2003) and/or a coronagraphic optical im- 
ager such as Eclipse (Trauger et al. 2000,2001) could obtain 
low-resolution spectra. 

To support the above efforts and planning for future programs 
of direct detection of extraoslar giant planets and to provide 
the theoretical context for the general analysis of the spectra 
and photometry of irradiated and isolated EGPs, our group has 
embarked upon a series of papers of EGP spectra, evolution, 
chemistry, transits, orbital phase functions, and light curves. 
The most recent paper in this series (SBH) explored generic 
features of irradiated EGP spectra as a function of orbital dis- 
tance, cloud properties, and composition class (Sudarsky, Bur- 
rows, and Pinto 2000). Our technical approach vis a vis radia- 
tive transfer, molecular abundance determinations, and cloud 
modeling is described in detail in that paper, to which the inter- 
ested reader is referred. However, SBH did not explore the di- 
agnostics of planetary mass and age. In addition, in their study 
of the orbital distance dependence of EGP spectra SBH did not 
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incorporate the effects of water or ammonia clouds in a fully 
consistent fashion. 

In the current paper, we allow our atmosphere code to de- 
termine cloud placement in a fully iterative, converged fashion. 
The result is a consistent determination of the dependence on 
distance of the spectra of EGPs irradiated by a G2V star, in- 
cluding the effects of the water and ammonia clouds that should 
form in their atmospheres. In §2, we summarize our numeri- 
cal approach. Then in §3 we present and describe our results 
for the dependence of the spectra of irradiated EGPs on orbital 
distance. In this paper, we emphasize the results for EGPs at 
wider separations (> 0.2 AU) and defer discussion of the cor- 
responding theory for close-in EGPs to a later paper^. This sec- 
tion treats the entire expected range of EGP emission/reflection 
spectra and behavior. In §4, we provide a representative se- 
quence of models that portray the dependence of irradiated EGP 
spectra on age, at a given mass and orbital distance. In §5, we 
present a representative sequence with EGP mass, at a given 
age and orbital distance. Section 6 is a digression into the ef- 
fect of stellar irradiation on companion brown dwarfs, charac- 
terized by much larger masses and more slowly decaying heat 
content. We study the signature in the optical of the reflection 
of stellar light from a companion brown dwarf. In fact, much 
of this paper is concerned with the signatures and diagnostics 
of the physical parameters of irradiated substellar-mass com- 
panions*. However, it is not feasible in one paper to explore 
all the possible combinations of planetary mass, age, compo- 
sition, orbital semi-major axis, eccentricity, and orbital phase 
with all stellar types. Hence, to maintain a reasonable focus, 
we restrict our discussions to G2V primaries, zero eccentricity 
orbits, and solar metallicity. We narrow our scope further by 
plotting only phase-averaged spectra (as in SBH) at zero orbital 
inclination. Papers on EGP orbital phase functions, albedos, 
and light curves are to foUow this one (e.g., Sudarsky, Burrows, 
and Hubeny 2004). These papers address the dependence on 
phase and Keplerian parameters. Finally, in §7 we present pre- 
dicted phase-averaged planet/star flux ratios for several known 
EGPs at wide separations and §8 reprises the essential conclu- 
sions of the paper. 

2. NUMERICAL ISSUES 

To calculate radiative/convective equilibrium atmospheres 
and spectra we use a specific variant of the computer pro- 
gram TLUSTY (Hubeny 1988; Hubeny & Lanz 1995). This 
variant involves the hybrid Complete Linearization/ Accelerated 
Lambda Iteration (CL/ALI) method, although in the present 
runs we use a full ALI mode which leads to an essential saving 
of computer time without slowing down the iteration process 
significantly. In addition, we employ TLUSTY's Discontinuous 
Finite Element (DEE) version (Castor, Dykema, & Klein 1992). 
The DEE technique, being first order, is optimal for handling ir- 
radiated atmospheres (SBH). All models are converged to one 
part in 10^. Stellar spectra from Kurucz (1994) are used for the 
incoming fluxes at the outer boundaries. The inner boundary 
condition is the interior flux and this (indexed by Tgff) is taken 
from the evolutionary models of Burrows et al. (1997) for the 
given mass and age, unless otherwise indicated. This approxi- 
mate procedure works well for wider- separation EGPs, but not 
as well for the closer-in EGPs (< 0. 15 AU). We use a standard 

' Note, however, that the spectra of close-in EGPs has been addressed in SBH, as 

and Goukenleuque et al. 2000. 

* We refer to Substellar-Mass Objects as SMOs. 



mixing-length prescription to handle convection and a mixing 
length of one pressure scale height. Since the atmosphere code 
is planar, we use the redistribution technique descibed in SBH, 
i.e. we weight the incident flux by 1/2 to account for the aver- 
age inclination of the planetary surface to the line of sight to the 
primary. As described in SBH, the planetary spectra we present 
here are phase/time-averaged over the orbit. Care is taken with 
this procedure to ensure that energy is conserved and that en- 
ergy in (from the star and the planetary interior) equals energy 
out. To account for the anisotropy of single scattering off of 
cloud particles, we calculate using Mie theory the average of 
the cosine of the scattering angle, and reduce the Mie-theory- 
derived total scattering cross section by one minus this average. 
With this procedure, we are substituting the "transport cross 
section" for the total scattering cross section. This approach 
has been shown to mimic the effect of asymmetric scattering 
quite well (Sudarsky, Burrows, and Pinto 2000). 

For molecular and atomic compositions, we use an updated 
version of the chemical code of Burrows and Sharp (1999), 
which includes a prescription to account for the rainout of con- 
densed species in a gravitational field and new thermochem- 
ical data. The derived molecular abundances are very simi- 
lar to those obtained by Lodders (1999) and Lodders and Fe- 
gley (2002). Our solar metalhcity is defined as the elemental 
abundance pattern found in Anders and Grevesse (1989). For 
molecular and atomic opacities, we have developed an exten- 
sive database, described in part in Burrows et al. (2001) and 
SBH. 

The treatment of H2O and NH3 clouds is done in a man- 
ner consistent with their respective condensation curves and the 
cloud base is put at the intersection of the corresponding con- 
densation curve at solar metallicity with the object's tempera- 
ture/pressure (T/P) profile. In each iteration of the global ALI 
scheme, we find a position of the cloud base as the intersec- 
tion of the corresponding condensation curve with the current 
T /P profile. The scale height of a cloud is assumed to be equal 
to one pressure scale height. We use Mie theory for the ab- 
sorptive and scattering opacities of particles whose model size 
is determined by the theory of Cooper et al. (2003). In the 
subsequent iteration of the ALI scheme, we employ this cloud 
opacity and scattering self-consistently in the radiative trans- 
fer equation and the energy balance equation. In the next ALI 
iteration we again recalculate the position of the cloud base, 
and the whole process is repeated until the cloud position is 
fully stabilized. We note that in the initial stages of the global 
iteration process the cloud position may vary significantly; in 
some case clouds appear (disappear) after several iterations of 
the cloudless (cloudy) atmosphere. This procedure ensures that 
the cloud position is self-consistent with the overall model at- 
mosphere. This also means that the predicted cloud properties, 
in particular optical depths and base pressures, vary in a physi- 
cal and consistent way with the orbital distance, mass, and age. 

3. ORBITAL DISTANCE DEPENDENCE OF EGP SPECTRA 
FROM 0.2 TO 15 AU 

Figure 1 shows the T/P profiles for the distance sequence 
from 0.2 to 15 AU for a 1-Mj EGP irradiated by a G2V star. For 
specificity, we have assumed a radius of 1 Rj and an internal flux 
Teff of 100 K for the entire family. According to the models of 

I as in Seager and Sasselov (1998,2000), Seager, Whitney, and Sasselov (2000), 
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Burrows et al. (1997), this corresponds roughly to an age of 5 
Gyr, but after ~0. 1 Gyr, the radius and gravity of the EGP vary 
little. The orbits are taken to be circular, so there is no assumed 
orbital phase dependence ^. The intercepts with the dashed lines 
identified by either { NH3 } or { H2O } denote the positions where 
the corresponding clouds form. Due to the cold trap effect and 
depletion due to rainout (Burrows and Sharp 1999), the higher- 
pressure intercept is taken to be at the base of the cloud. The 
spectral/atmospheric models include the effects of these clouds 
in a consistent way. Table 1 gives the modal particle sizes in 
microns that we derive using the theory of Cooper et al. (2003) 
when a cloud of either water or ammonia (or both) appears. For 
the ammonia clouds that form in this orbital distance sequence 
(at >6 AU), the modal particle sizes we find hover near 50-60 
/xm . The corresponding particle sizes in the water clouds are 
near 110 fim . These particles are larger than for more massive 
SMOs. Furthermore, we assume that the particle size does not 
vary with altitude and that the particle size distribution (given 
a modal radius) is that of Deirmendjian (1964,1969). Clearly, 
a major ambiguity in EGP modeling is cloud physics (SBH). 
We have settled on the Cooper et al. (2003) theory to provide a 
consistent framework. 

When the chemistry indicates that both cloud types are 
present, we include them both in the atmospheric/spectral cal- 
culation. As Fig. 1 shows, the ammonia cloud is always above 
the water cloud. Water clouds form around a G2V star exte- 
rior to a distance near 1 .5 AU, whereas ammonia clouds form 
around such a star exterior to a distance near 4.5 AU. Note that 
Jupiter itself is at the distance from the Sun of ~5.2 AU. For the 
closer-in objects, the T/P profile manifests an inflection. This 
inflection is a consequence of the dominance of external radia- 
tion over the internal heat flux. For longer ages and low masses, 
the orbital distance at which one must place the SMO/EGP to 
erase this inflection is large. For this model set, that distance is 
~2 AU. For larger masses, dimmer primaries, and shorter ages, 
that distance decreases. In addition, the strong irradiation that 
produces the inflection in the T/P profile also forces the radia- 
tive/convective boundary to recede to higher pressures. For a 
1-MjEGP at an orbital distance of 0.05 AU (not shown) around 
a G2V star, this pressure can be greater than 1000 bars! Such 
is the case for HD 209458b and OGLE-TR56b (Burrows, Su- 
darsky, and Hubbard 2003; Fortney et al. 2003). 

For comparison, and in anticipation of the discussion in §4 
concerning the age dependence of irradiated EGP spectra. Fig. 
2 portrays the evolution in the T/P profile of a 1-MjEGP at a 
distance of 4 AU around a G2V star. Table 2 gives the corre- 
sponding modal particle sizes in microns (Cooper et al. 2003) 
for the ice particles of the water clouds that appear in this evo- 
lutionary sequence, as well as the inner flux Teff, gravity, and 
planet radius as this 1-MjEGP evolves according to the the- 
ory of Burrows et al. (1997). The modal particle size is very 
roughly constant with age. The radius of the 1-MjEGP de- 
creases by ~15% from 0.1 to 5 Gyr. Not surprisingly, at a 
distance of 4 AU no inflection in the T/P profile is produced. 
Note that at this orbital distance and as early as ~50 Myr (not 
shown), water clouds form in Jovian-mass objects. Note also 
that at 4 AU, even after 5 Gyr ammonia clouds have not yet 
formed in the atmosphere of an irradiated 1 -Mj EGP. This is not 
true for a similar object in isolation (Burrows, Sudarsky, and 
Lunine 2003). 

Figure 3 depicts the planet- to- star flux ratios from 0.5 /umto 



30 ^mfor the orbital distance study associated with the T/P 
profiles shown in Fig. 1. In the optical, the flux ratios vary 
between 10"^ and 10"^*^. In the near infrared, this ratio varies 
widely from ~10~^ to 10"^*. However, in the mid-infrared be- 
yond 10 yum, the flux ratio varies more narrowly from 10"^ to 
10"^. Hence, it makes a difference in what wavelength region 
one conducts a search for direct planetary light. 

Gaseous water absorption features (for all orbits) and 
methane absorption features (for the outer orbits) sculpt the 
spectra. The reflected component due to Rayleigh scattering 
and clouds (when present) is most manifest in the optical and 
the emission component (similar to the spectrum of an isolated 
low-gravity brown dwarf) takes over at longer wavelengths. For 
the EGPs interior to ^1.0 AU, the fluxes longward of ^0.8 
/xmare primarily due to thermal emission, not reflection. These 
atmospheres do not contain condensates and are heated effi- 
ciently by stellar irradiation. As a result, the Z (~1.0 /itm), 
J (^1.2 fim), H (^1.6 yum), and K (^2.2 /im) band fluxes 
are larger by up to several orders of magnitude than those of 
the more distant EGPs. Generally, clouds increase a planet's 
flux in the optical, while decreasing it in the /, K, L' (~3.5 
/im), and M (~5.0 /xm) bands. The transition between the 
reflection and emission components moves to longer wave- 
lengths with increasing distance, and is around 0.8-1.0 /xmat 
0.2 AU and ~3.0 /xmat 15 AU. However, since the irradiation 
and atmospheric structure and spectra are being calculated self- 
consistently, emission and reflection components are in fact in- 
extricibly intertwined and it is not conceptually correct to sep- 
arate them. 

At large distances exterior to ~3.5 AU, the flux longward 
of ^15 /im manifests undulations due to pressure-induced ab- 
sorption by H2. Importantly, there is always a significant bump 
around the M band at 4-5 /im . The peak of this bump shifts 
from ~4 /imto ~5 /xm with increasing orbital distance. Though 
muted by the presence of clouds, it is always a prominent fea- 
ture of irradiated EGPs, as it is in T dwarfs and the Jovian plan- 
ets of our solar system. Curiously, but not unexpectedly, as 
Fig. 3 indicates, the planet/star flux ratio is most favorable in 
the mid-infrared. This fact should be of some interest to those 
planning TPF or successor missions to the Spitzer Space Tele- 
scope. 

Though the major trend is a monotonic decrease in a planet's 
flux with increasing orbital distance, the 0.2-AU model at- 
mosphere is hot enough that the sodium and potassium reso- 
nance absorption lines appear and surpress the flux around Na- 
D (0.589 /im) and die related K I doublet at 0.77 /im. The 
result is a lower integrated visible flux that is comparable to 
that of the otherwise dinnmer 0.5-AU model. Figures 4 and 5 
focus in on the 0.5 /imto 2.0 /imregion and allow one to dis- 
tinguish one model from another at shorter wavelengths more 
easily than is possible in the panoramic Fig. 3. These figures al- 
low us to see that for greater orbital distances, the atmospheric 
temperatures are too low for the alkali metals to appear, but the 
methane features near 0.62 /im, 0.74 /im, 0.81 /im, and 0.89 
/imcome into their own. Broad water bands around 0.94 /im, 
1.15 /im, 1.5 /imand 1.85 /imthat help to define the Z, J, and 
H bands are always in evidence, particularly for the 0.2 and 
1.0 AU models that don't contain water clouds. For greater dis- 
tances, the presence of water clouds slightly mutes the variation 
with wavelength in the planetary spectra. Hence, smoothed wa- 
ter features and methane bands predominate beyond ~1.5 AU. 



^ Note that with a significant eccentricity, this assumption is not valid. 
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4. AGE DEPENDENCE FROM 100 MYR TO 5 GYR OF THE 

SPECTRUM OF A 1-MjEGP AT A GIVEN DISTANCE 

Figure 6 presents the planet-to-star flux ratios from 0.5 /imto 
6.0 /umfor a 1-MjEGP orbiting a G2V star at 4 AU as a func- 
tion of age. Figure 2 depicts the corresponding T /P profiles, 
along with the NH3 and H2O condensation lines at solar metal- 
licity. The theory of Burrows et al. (1997) is used to obtain 
an approximate mapping between mass, age, internal Tgff, and 
gravity^. The depicted ages are 0.1, 0.3, 1, 3, and 5 Gyr. These 
planet parameters are chosen merely to represent the system- 
atics with age; different EGP masses and orbital distances will 
yield quantitatively different spectra. Table 2 shows that for this 
suite of models the internal flux Tgff varies from 290 K to 103 
K, the surface gravity varies from 1695 cm s"'^ to 2325 cm s"'^, 
and the planet radius varies from 1 . 17 /?j to 1 .0 /?j . 

As is clear from Fig. 6, younger EGPs with higher inner 
boundary TgffS (Table 2) have much higher fluxes in the Z, J, 
H, K, and M bands. However, the older objects, having cooled 
more, have lower internal luminosities. This results in lower 
fluxes in those same bands by as much as two orders of mag- 
nitude. As Figure 2 shows, the older EGPs have progressively 
deeper water clouds. What is not obvious from Figure 2 is that 
these clouds are also thicker. This results in a very slightly in- 
creasing reflected optical flux with increasing age that accom- 
panies the reverse trend in the near infrared. Hence, the fluxes 
in the optical shortward of ^1.0 /imare only weak functions of 
age, while the fluxes in the Z, J, H, K, and M bands at early 
ages are strong functions of age. At later ages, the formation 
of water clouds moderates the age dependence of the Z, J, H, 
and K band fluxes. In fact, there can be slight increases in flux 
in the Z, J, and H bands with increasing age. However, the M 
band flux continues to be diagnostic of age, monotonically de- 
creasing by almost two orders of magnitude from 0.1 to 5 Gyr. 
Hence, the best diagnostics of age are in the near infrared, not 
the optical. 

5. IRRADIATED EGP SPECTRA AS A FUNCTION OF MASS 

FROM 0.5MjTO 8MjAT A GIVEN AGE AND ORBITAL 

DISTANCE 

Figure 7 portrays planet-to-star flux ratios from 0.4 ;umto 6.0 
/zmfor a 5-Gyr EGP orbiting a G2V star at 4 AU, as a function 
of EGP mass. The masses represented are 0.5, 1, 2, 4, 6, and 8 
Mj . An inner flux boundary condition Tgff from the evolution- 
ary calculations of Burrows et al. (1997) has been employed 
and is given in Table 3. T^s varies from 82 K to 251 K, the 
surface gravity varies from 1290 cm s"-^ to 17800 cm s"-^, and 
the radii vary from 0.95 Rj to 1 .04 Rj . Note that these radii 
peak in the middle of the sequence near 4 Mj . Due to the 5- 
Gyr age assumed, all these models have water clouds and the 
derived modal particle sizes decrease monotonically with in- 
creasing mass from 146 /zmat 0.5 Mj to 39 /imat 8 Mj . In gen- 
eral, the larger the EGP mass, the higher in the atmosphere the 
clouds form, but cloud position is not fully monotonic at the 
low-mass end of the sequence. 

For higher mass, at a given age an EGP's inner Tgff and in- 
ternal luminosity are higher. This results in higher fluxes in the 
Z, J, H, K, and M bands for higher masses and is similar to 
the trend seen in §4 with decreasing age. However, larger-mass 
EGPs also have higher surface gravities, which result in wa- 
ter clouds with lower column depths, and, hence, lower optical 

^ We note that the star evolves as well, but for clarity we have neglected this effect. 
' For comparison, some of the closest EGPs are also shown at the bottom of Table 4. 



depths, despite the contrary trend of modal particle size (Table 
3). The upshot is that higher-mass EGPs have sUghtly smaller 
planet/star flux ratios in the optical. This optical component is 
due to reflection off of cloudy atmospheres with roughly similar 
compositions. Hence, there is an anti-correlation between flux 
levels in the visible and near-infrared that might be diagnostic 
of planet mass for Jupiter-aged EGPs. 

6. IRRADIATED BROWN DWARFS 

Depicted in Fig. 8 are theoretical spectra of a 30-Mj brown 
dwarf at ages of 1 and 5 Gyr in orbit around a G2V star from 
0.4 /imto 1.5 iJ,m. The results are for orbital distances from 
5 AU to 40 AU and a distance to the Earth of 10 parsecs and 
include irradiation effects. As before, the theory of Burrows et 
al. (1997) is used to determine Teff and gravity for this mass 
and these ages. In Fig. 8, the prominence of the Na-D (0.589 
/xm) and K I (0.77 /xm) features at shorter wavelengths is clear 
and is canonical for brown dwarfs (Burrows, Marley, and Sharp 
2000). 

Unlike for the lower mass EGPs discussed in §3, §4, and §5, 
the internal luminosity of such a relatively massive SMO dom- 
inates its energy budget. As a consequence, the brown dwarf's 
spectrum longward of 0.9 /imis unaffected by irradiation. How- 
ever, the reflected component in the optical, particularly for the 
older brown dwarf with lower internal heat content and lumi- 
nosity, is a function of distance. As Fig. 8 demonstrates, the 
optical flux from a cool brown dwarf can be elevated shortward 
of 0.7 /imby as much as an order of magnitude in the V band. 
In particular, the shape of the Na-D feature can be significantly 
altered. Since there are no clouds in the atmospheres of these 
brown dwarf models, Rayleigh scattering off of H2, He, and 
H2O accounts for this reflection. Our prediction is that the opti- 
cal spectra of brown dwarfs that are close companions to K, G, 
or F stars will be modified by irradiation. Note that the T dwarf 
Gliese 229B is at a projected distance of ^^40 AU, but that its 
primary is an M4V star. Such a star is too dim to so radically 
alter a brown dwarf's optical flux. For an L dwarf companion, 
the presence of siUcate clouds in its atmosphere will reflect a 
primary's light in the blue and UV. Someday, such a reflected 
component might be detectable. 

7. PREDICTED PHASE- AVERAGED SPECTRA FOR KNOWN 

EGPS AT WIDE ANGULAR SEPARATIONS 

Table 4 Usts many of the known EGPs that, due to a pro- 
pitious combination of semi-major axis and distance from the 
Earth, are at wide angular separations from their parent stars^. 
Ordered by decreasing separation (defined as the ratio of semi- 
major axis to distance). Table 4 also lists the stellar type of the 
primary, semi-major axis, Hipparcos distance, orbital period, 
nip sin(;), and orbital eccentricity. This family of known EGP 
comprises some of the prime candidates for direct detection of 
planetary Ught using interferometric, adaptive-optics, or coro- 
nagraphic techniques (§1). Note that the separations quoted 
in Table 4 don't take into account the projection of the or- 
bit or variations due to non-zero eccentricities (which can be 
large). In particular, variations due to the significant excursions 
in planet-star distance that attend large eccentricities can result 
in large changes in irradiation regimes. In turn, this can result in 
large variations in planet spectrum. Interestingly, it is possible 
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for an EGP atmosphere to cycle between states with and with- 
out clouds, with the concomitant large changes in flux ratios 
and spectral signatures with orbital phase (Sudarsky, Burrows, 
and Hubeny 2004). Figure 3 gives some idea of the range of 
spectral variation possible for highly-eccentric EGPs. 

Figures 3, 6, and 7 provide a broad-brush view of generic 
EGP spectra for a range of orbital distances, ages, and masses. 
With Figs. 9, 10, and 11, we provide phase-averaged predic- 
tions/calculations for a specific subset of the known EGPs listed 
in Table 4. This subset includes HD 39091b, 7 Cephei b, HD 
70642b, V And d, Ghese 777A b, HD 216437b, HD 147513b, 
55 Cancri d, 47 UMa b, 47 UMa c, 14 Her b, and e Eri b. Table 5 
itemizes these EGPs, along with their derived TeffS and surface 
gravities. In addition. Table 5 lists the theoretical modal particle 
sizes of the water droplets that form in their atmospheres. For 
definiteness, the planetary masses are set equal to the measured 
Mp sin(j) and, as before, the evolutionary theory of Burrows et 
al. (1997) is used to estimate the corresponding surface gravi- 
ties and inner boundary TgffS. As described in SBH and §2, a 
stellar spectrum from Kurucz (1994) for the primary stellar type 
given for each EGP in Table 5 is used for the outer irradiation 
boundary condition in the self-consistent atmosphere/spectrum 
calculation. As a comparison of Figs. 3, 4, and 5 with Figs. 9, 
10, and 11 demonstrates, the major dependence of the planet- 
to-star flux ratio is with orbital distance. Note that for aU of 
the EGPs listed in Table 5, the assumption that the planet-star 
distance is fixed at the measured semi-major axis does not re- 
sult in the formation of ammonia clouds. However, such clouds 
should form near apastron for the known EGPs with large semi- 
major axes and high eccentricities, such as e Eri b and 55 Cnc d. 
The appearance and disappearance of such clouds with orbital 
phase would be exciting signatures to detect. 

8. SUMMARY AND CONCLUSIONS 

In this paper, we have calculated theoretical phase-averaged 
planet/star flux ratios in the optical, near-infrared, and mid- 
infrared for wide-separation irradiated EGPs as a function of 
orbital distance, mass, and age. We have also predicted the cor- 
responding quantities for 12 specific known EGPs, given their 
corresponding primary spectra, average orbital distances, ap- 
proximate masses, and approximate ages. Hence, we have ex- 
plored various physical diagnostics that can inform the direct 
EGP detection programs now being planned or proposed. In 
the optical, the flux ratios for distances from 0.2 AU to 15 AU 
and masses from 0.5 Mjto 8 Mjcan vary from slightly above 
10"** to ^10"'". In the near infrared around 2-4 /im, the flux 
ratio ranges more widely, spanning values from 10"'* to as low 
as 10"*^. At ~ 5 ^m, the planet/star flux ratio ranges approx- 
imately four orders of magnitude and can be as high as lO"'*. 
The M band should be a useful region to explore and M-band 
fluxes are sensitive to the presence of clouds. For all models, 
the mid-infrared from 10 /umto 30 /umis always encouragingly 
high. In fact, for closer separations (0.05 AU-0.1 AU), not the 
subject of this paper, the flux ratio at ^20 //mean be ^10"^. 

Depending upon orbital distance, age, and mass, spectral fea- 
tures due to methane, water, and the alkali metals are promi- 
nent. Furthermore, there is a slight anti-correlation in the ef- 
fects of clouds in the optical and infrared, with the optical fluxes 



increasing and the infrared fluxes decreasing with increasing 
cloud depth. For young and massive irradiated EGPs, there are 
prominent peaks in the Z, J, H, K, and M bands. Though the 
optical flux is not a very sensitive function of age, there is a use- 
ful age dependence of the fluxes in these bands. Furthermore, 
there is an anti-correlation with increasing mass at a given age 
and orbital distance between the change in flux in the optical 
and in the near-IR bands, with the optical fluxes decreasing and 
the Z, /, H, K, and M band fluxes increasing with increasing 
mass. 

The spectra of more massive SMOs (brown dwarfs) long- 
ward of ~0.9 jumare not significantly affected by stellar irra- 
diation. Their internal heat content and interior fluxes are too 
large. However, brown dwarf fluxes from 0.4 /xmto 0.65 /xm, 
can be enhanced by Rayleigh reflection by as much as a factor 
of 10. This is particularly true of old or low-mass brown dwarfs 
and is a predictable function of distance. Moreover, irradiation 
can alter the profile shape of the Na-D feature significantly. 

We are not able to calculate EGP spectra for all possible com- 
binations of mass, age, composition, orbital distance, eccentric- 
ity, orbital phase, Keplerian element, and primary spectral type. 
This fact is what motivates the more modest synoptic view we 
have provided in this paper However, in the process of devel- 
oping the tools for this study, we have estabhshed the capabil- 
ity to calculate EGP spectra for any combination of these pa- 
rameters. In particular, Sudarsky, Burrows, and Hubeny (2004) 
address the orbital phase and eccentricity dependences of irra- 
diated EGPs spectra. 

The remote sensing of the atmospheres of EGPs will be chal- 
lenging, but the detection and characterization of the direct Ught 
from a planet outside our solar system will be an important 
milestone in both astronomy and planetary science. One instru- 
ment proposed to meet this challenge is the space-based corona- 
graphic imager Eclipse (Trauger et al. 2000,2001). The Eclipse 
instrument team is predicting a contrast capability of ~10"^ for 
an inner working angle of 0.3" or 0.46" in the V/R or Z bands, 
respectively. As Figs. 3 through 7 indicate, with such a capa- 
biUty irradiated EGPs could be detected and analyzed. Many 
could even be discovered. However, whether such sensitivity 
is achievable remains to be demonstrated. Be that as it may, 
further advancement in our understanding of extrasolar plan- 
ets is contingent upon technical advances that would enable the 
direct measurement of the dynamicaUy dominant and brighter 
components of extrasolar planetary systems, the EGPs. 
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Table 1 



Model Data for Distance Sequence 



a (AU) 


condensale(s) 


I'o (pm) 


0.2 


none 




0.5 


none 




1 


none 




2 


H2O 


107 


4 


H2O 


109 


6 


NH3, H2O 


56, 110 


8 


NH3, H2O 


56, 111 


10 


NH3, H2O 


57, 111 


15 


NH3, H2O 


62, 111 



Table 2 



Model Data for Age Sequence 



Age (Gyr) 


Teff(K) 


logw 8 (cm s ^) 


R/Rj 


ro (/um) 


0.1 


290 


3.23 


1.17 


123 


0.3 


219 


3.27 


1.11 


119 


1 


159 


3.31 


1.06 


123 


3 


118 


3.35 


1.02 


111 


5 


103 


3.37 


1.00 


109 






Table 3 








Model Data for Mass Sequence 


M/M] 


Teff(K) 


logio 8 (cm s 2) 


R/Rj 


ro (pm) 


0.5 


82 


3.11 


0.95 


146 


1 


103 


3.37 


1.00 


109 


2 


134 


3.64 


1.03 


80 


4 


177 


3.93 


1.04 


57 


6 


216 


4.12 


1.03 


46 


8 


251 


4.25 


1.02 


39 



Table 4 



Interesting EGPs Listed by Angular Separation 



EGP 


separation (") 


star 


a(AU) 


d(pc) 


P 


Msin; (Mj ) 


e 


e Eri b 


1.0 


K2V 


3.3 


3.2 


6.85 yrs. 


0.86 


0.61 


55 Cnc d 


0.44 


G8V 


5.9 


13.4 


14.7 


4.05 


0.16 


47 UMa c 


0.28 


GOV 


3.73 


13.3 


7.10 


0.76 


0.1 


Gl 777A b 


0.23 


G6V 


3.65 


15.9 


7.15 


1.15 


~0 


V And d 


0.19 


F8V 


2.50 


13.5 


3.47 


4.61 


0.41 


HD 39091b 


0.16 


GlIV 


3.34 


20.6 


5.70 


10.3 


0.62 


47 UMab 


0.16 


GOV 


2.09 


13.3 


2.98 


2.54 


0.06 


7 Cephei b 


0.15 


K2V 


1.8 


11.8 


2.5 


1.25 


~0 


HD 160691c 


0.15 


G3IV-V 


2.3 


15.3 


3.56 


~1 


~0.8 


14 Her b 


0.15 


KOV 


2.5 


17 


4.51 


3.3 


0.33 


HD 33636b 


0.12 


GOV 


3.56 


28.7 


4.43 


7.71 


0.41 


HD 10647b 


0.12 


F9V 


2.10 


17.3 


2.89 


1.17 


0.32 


HD 70642b 


0.11 


G5IV-V 


3.3 


29 


4.79 


2.0 


0.10 


HD 216437b 


0.10 


G4V 


2.7 


26.5 


3.54 


2.1 


0.34 


HD 147513b 


0.098 


G3V 


1.26 


12.9 


1.48 


1.0 


0.52 


HD 160691b 


0.097 


G3IV-V 


1.48 


15.3 


1.74 


1.7 


0.31 


HD 168443c 


0.087 


G5V 


2.87 


33 


4.76 


17.1 


0.23 


HD 50554b 


0.077 


F8V 


2.38 


31.03 


3.50 


4.9 


0.42 


HD 106252b 


0.070 


GOV 


2.61 


37.44 


4.11 


6.81 


0.54 


HD 10697b 


0.067 


G5IV 


2.0 


30 


2.99 


6.59 


0.12 


V And c 


0.061 


F8V 


0.83 


13.5 


241 days 


2.11 


0.18 


GJ 876b 


0.045 


M4V 


0.21 


4.72 


61.0 


1.89 


0.1 


GJ 876c 


0.028 


M4V 


0.13 


4.72 


30.1 


0.56 


0.27 


HD 114762b 


0.013 


F9V 


0.35 


28 


84.0 


11.0 


0.34 


55 Cnc b 


8.2 X 10"^ 


G8V 


0.12 


13.4 


14.7 


0.84 


0.02 


V And b 


4.4 X 10-3 


F8V 


0.059 


13.5 


4.62 


0.71 


0.034 


51 Peg b 


3.4 X 10-3 


G2V 


0.05 


14.7 


4.23 


0.44 


0.01 


T Boo b 


3.3 X 10-3 


F7V 


0.05 


15 


3.31 


4.09 


~0 


HD 209458b 


9.6 X 10-* 


GOV 


0.045 


47 


3.52 


0.69 


~0 


HD 83443b 


8.7 X 10-* 


KOV 


0.038 


43.5 


2.99 


0.35 


0.08 
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Table 5 

Model Data for Specific EGPs 



EGP 


star 


a(AU) 


d(pc) 


Msinj 


Teff(K) 


logm g (cgs) 


ro (/xm) 


e Eri b 


K2V 


3.3 


3.2 


0.86 


170 


3.48 


51 


55 Cnc d 


G8V 


5.9 


13.4 


4.05 


200 


4.30 


51 


47 UMa c 


GOV 


3.73 


13.3 


0.76 


95 


3.48 


96 


G1777Ab 


G6V 


3.65 


15.9 


1.15 


100 


3.48 


97 


V And d 


F8V 


2.50 


13.5 


4.61 


250 


4.30 


51 


HD 39091b 


GlIV 


3.34 


20.6 


10.3 


300 


4.48 




47 UMa b 


GOV 


2.09 


13.3 


2.54 


150 


3.78 


67 


7 Cephei b 


K2V 


1.8 


11.8 


1.25 


150 


3.60 


83 


14 Her b 


KOV 


2.5 


17 


3.3 


150 


3.90 


59 


HD 70642b 


G5V 


3.3 


29 


2.0 


130 


3.64 


80 


HD 216437b 


G4V 


2.7 


26.5 


2.1 


150 


3.70 


72 


HD 147513b 


G3V 


1.26 


12.9 


1.0 


100 


3.37 


96 



References. — 1) From Sudarsky, Burrows, & Hubeny (2003), where a modal particle size of 5 /xm was 
assumed. 
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1500 



T [K] 



Fig. 1 . — Profiles of atmospheric temperature (in Kelvin) versus the logarithm base ten of the pressure (in bars) for a family of in'adiated 1-Mj EGPs around a 
G2V star as a function of orbital distance. Note that the pressure is decreasing along the ordinate, which thereby resembles altitude. The orbits ai'e assumed to be 
circular, the planets are assumed to have a radius of 1 Sj , the effective temperature of the inner boundary flux is set equal to 100 K, and the orbital separations vary 
from 0.2 AU to 15 AU. The intercepts with the dashed lines identified with either {NHj} or {H2O} denote the positions where the corresponding clouds form. See 
text for a discussion. 
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1 1 Mj planet orbiting a G2V star at 4 AU 
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T [K] 



Fig. 2. — Similar to Fig. 1. this figure depicts atmosplieric profiles of temperature (in Kelvin) versus the logarithm base ten of the pressure (in bars) for a family 
of irradiated i-Mj EGPs around a G2V star, but as a function of age at a given orbital distance of 4.0 AU. The ages vary from 0. 1 to 5.0 Gyr. Note that the pressure 
is decreasing along the ordinate, which thereby resembles altitude. The cloud condensation curves for both NH3 (moot for this sequence) and HtO are given as the 
dashed lines and the spectral/atmospheric models include the effects of the water clouds in a consistent way. See text in §3 and §4 for relevant details and discussion. 
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Wavelength [microns] 



Fig. 3. — Planet to star flux ratios versus wavelength (in microns) from 0.5 fimto 30 fimfor a 1-Mj EGP with an age of 5 Gyr orbiting a G2V main sequence star 
similar to the Sun. This figure portrays ratio spectra as a function of orbital diatance from 0.2 AU to 15 AU. Zero eccentricity is assumed and the planet spectra have 
been phase-averaged as described in Sudarsky, Burrows, and Hubeny (2003). The associated T/P profiles are given in Fig. 1 and Table 1 lists the modal radii for 
the particles in the water and ammonia clouds. Note that the planet/star flux ratio is most favorable in the mid-infrared. See text for further discussion. 
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Wavelength [microns] 



Fig. 4. — The same as Fig. 3, but highlighting the shorter wavelengths and for a subset of distances (0.2, 1,4, 10 AU). This figure provides a clearer picture of the 
features shortward of 2.0 /jmfor each of the represented models. For the 0.2 AU model, the temperatures of the atmosphere are high enough for the Na-D doublet 
around 0.589 fimand the K I doublet near 0.77 fj,mto be visible. These features are even more prominent for closer-in EGPs (Sudarsky, Burrows, and Hubeny 2003). 
At greater orbital distances, the atmospheric temperatures are too low for the alkali metals to appear, but the methane features near 0.62 fim , 0.74 fim , 0.8 1 fj,m , and 
0.89 /jmcome into their own. Water bands around 0.94 fim, 1.15 fim, 1.5 fimand 1.85 fimthaX help to define the Z, J, and H bands are always of importance. For 
greater distances, the presence of water clouds can smooth the variations in the planetary spectra that would otherwise be large due to the strong absoiption features 
of gaseous water vapor. See text for discussion. 
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Fig. 5. — The same as Fig. 3. but, as in Fig. 4, liigliligliting the sliorter wavelengths. A different subset of distances (0.5, 2, 6, 15 AU) is shown. See the text and 
the figure caption for Fig. 4 for details. Figures 4 and 5 allow one to distinguish more easily than is possible in the panoramic Fig. 3 one model from another. 
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Fig. 6. — The planet-to-star flux ratio from 0.5 /jmto 6.0 fimfoT a l-Mj EGP orbiting a G2V star at 4 AU as a function of age. The ages are 0.1, 0.3, 1, 3, and 5 
Gyr. An inner flux boundary condition Tefffrom the evolutionary calculations of Burrows et al. (1997) has been employed. The effect of clouds is handled in the 
radiative transfer calculation in a completely consistent fashion. See Table 2, Figure 2, and text for details and discussion. 
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Fig. 7. — Similar to Fig. 6, but tlie planet-to-star flux ratio from 0.4 ^mto 6.0 ^mfor a 5-Gyr EGP orbiting a G2V star at 4 AU, as a function of EGP mass. The 
masses represented are 0.5, 1, 2, 4, 6, and 8 Mj . See Table 3 and text for a discussion. 
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A [^im] 



Fig. 8. — Depicted are theoretical spectra of a 30-Mj brown dwarf at ages of 1 and 5 Gyr in orbit around a G2V star that is irradiating it. The logarithm base ten 
of the flux in milliJanskys at 10 parsecs versus wavelength in microns from 0.4 fimto 1.5 ^mis given. The theory of Burrows et al. (1997) was used to determine 
Teff and gravity for these ages and mass. The results are shown for different orbital distances from 5 to 40 AU. See text for a discussion. 
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Fig. 9. — Similar to Figs. 6 and 7, but for four specific known EGPs listed in Table 5. These are HD 39091b, 7 Cephei b, HD 70642b, and v And d. The 
wavelengths range from 0.5 /imto 6.0 fim. The planet fluxes are phase-averaged and the effect of the known eccentricities is ignored. The orbital distances are 
assumed to be equal to the measured semi-major axes and the planets' masses are set equal to the measured values of nip s'm(i). For all of these objects water clouds 
are formed and the modal particle sizes are given in Table 5. Table 5 also shows the values of the gravity and the T^ff that describes the inner flux, both derived using 
the theory of Burrows et al. (1997). Compare this figure with Figs. 3, 6, and 7 and see text for details. 
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Wavelength [microns] 



Fig. 10.— Same as Fig. 9, but for Gliese 777A b, HD 216437b, HD 147513b, and 55 Cancri d. Refer to Fig. 9, Table 5, and the text for further details and 
discussion. 
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Wavelength [microns] 



Fig. 1 1 . — Same as Figs. 9 and 10, but for 47 UMa b, 47 UMa c, 14 Her b, and e Eridani b. Refer to Fig. 9, Table 5, and the text for further details and discussion. 



